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Band structures of Cu2ZnSnS4 and Cu2ZnSnSe4 from many-body methods 
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We calculate the band structures of kesterite and stannite Cu2ZnSnS4 and CU2Z11S11SC-4, using a state-of- 
the-art self-consistent GW approach. Our accurate quasiparticle states allow to discuss: the dependence of 
the gap on the anion displacement; the key-role of the non-locality of the exchange-correlation potential to 
obtain good structural parameters; the reliability of less expensive hybrid functional and GGA+U approaches. 
In particular, we show that even if the band gap is correctly reproduced by hybrid functionals, the band- 
edge corrections are in disagreement with self-consistent GW results, which has decisive implications for the 
positioning of the defect levels in the band gap. 



Thin-film solar cells made of Cu(In,Ga)Se 2 (CIGS) ally 
cost reduction and high efficiency, and compete today as 
successors of the dominating silicon technology. Never- 
theless, there are concerns about their large scale pro- 
duction due to the increasing price of In and Ga. Qua- 
ternary chalcogenidcs Cu2ZnSn(S,Se)4 (CZTS) have re- 
cently been proposed as alternative absorbers. Their 
crystal structures and electronic properties are very sim- 
ilar to those of the parent CIGS, while their constituent 
elements are naturally abundant and non-toxic. The al- 
loys C^ZnSnSajSci-z have optimal gaps according to the 
Shockley-Queisser limit and their use as absorbers in thin 
film solar cells is getting established by a growing energy 
conversion efficiency (almost 10%^ for lab cells). How- 
ever, the understanding of the properties of the different 
phases of Cu2ZnSnS4 and Cu2ZnSnSe4is still rather su- 
perficial, and only few recent studies have addressed their 
structural,— ~— electronic,—^ and defect properties.— £ 

The zincblende-derived kesterite structure (14) of 
CZTS is recognized to be the most stable. The energy 
difference per atom with respect to the stannite struc- 
ture (/42m) is only of few mcV per atom^^ proving 
that kesterite and stannite phases can coexist in ex- 
perimental samples, and explaining the reported disor- 
dered structures^ Similarly to chalcopyrites CIGS, qua- 
ternary kesterite and stannite CZTS are obtained from 
the zincblende structure by replacing the Zn cations in 
such a way that each anion (Se or S) is coordinated by 
one Zn, one Sn and two Cu atoms. The existence of 
three distinct cations results in three different cation- 
anion bond lengths, which induce a displacement of the 
anion from its ideal zincblende site. That distortion is 
measured by the anion displacement parameters (u x , u y , 
u z ), i.e., the relative coordinates of the anion in the con- 
ventional body-centered tetragonal cell. The anion dis- 
placement is harder to measure than the lattice constants 
a and c, due to the inhomogencity of the samples. In- 
deed, the dispersion of data for u is known to be large 
in CIGS compounds ! 11 ^ 2 Concerning CZTS, fewer mea- 
surements can be found in literature and in most cases 
only a and c are reported^ Concerning the experimental 
band gaps,— ~— early studies suggested a value of about 
1.4-1.6 eV both for S and Se compounds. This was in 
disagreement with density functional theory (DFT) cal- 




FIG. 1. Band structures for kesterite Cu2ZnSnS4 calculated 
using: (left panel) DFT-LDA (dashed line) and scGW (con- 
tinuous line), (right panel) GGA+U (dashed line) and HSE06 
(continuous line). Top valence bands are set to zero. 



culations,— & using semi-local or hybrid functionals, which 
obtained systematically a smaller gap for Se compounds. 
More recent measurements^ delineate a gap of about 
leV for Cu2ZnSnSe4, explaining the previous overesti- 
mation with the presence of ZnSc in the sample. 

Note that the current workhorse ab initio theory is, 
in fact, DFT, in combination with the local density ap- 
proximation (LDA) or generalized gradient approxima- 
tions (GGA). However, this approach is totally inade- 
quate to study the electronic structure of materials where 
the band gap is controlled by the hybridization of the d 
states of a transition metal with p states, such as in CIGS 
or Cu2ZnSnS4. Moreover the anion-cation bonds can be 
poorly described, leading to anion displacements outside 
the experimental range, even if lattice constants are as 
usual well reproduced. In particular, for CuIn(S,Se)2, 
LDA and GGA yield structural internal parameters with 
an error of about 5% for it . 11 ! 12 ' 17 Such an error leads to 
an underestimate by more than 50% of the band gap in 
a calculation at the theoretical geometry ji 7 - 
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In this letter, we solve these problems by using state- 
of-the-art ab initio approaches that go beyond standard 
DFT to determine accurate quasiparticle band struc- 
tures of both kesterite and stannite Cu2ZnSnS4 and 
Cu2ZnSnSe4. These calculations are based on a re- 
stricted self-consistent (sc) GW scheme, which has the 
advantage of being independent of the starting point (i.e., 
the poor LDA Kohn-Sham states) at the price of a larger 
computational complexity. Such approach, that we will 
refer to as scGW, consists in performing a self-consistent 
GWcalculation within the Coulomb hole plus screened 
exchange (COHSEX) approximation^ followed by a per- 
turbative GW on top of it. This method has been applied 
to many transition-metal compounds, yielding excellent 
results for the band gaps and the quasiparticle band 
structure Standard LDA or GGA, and scGW cal- 

culations were performed using the code ABINIT. We 
included semicore states in the valence to build the norm- 
conserving pseudopotentials for Cu, Zn and Sn. We 
also used the code VASP— for Heyd-Scuseria-Ernzerhof 
(HSE06)22 hybrid functional and GGA+U calculations. 
Due to the similarity of CZTS materials with the CIGS 
family, the convergence parameters turned out to be the 
same reported in Rcfs. [l7l 

In Fig. [T] we display band structures for kesterite 
Cu2ZnSnS<4, obtained using different theoretical schemes 
at the experimental geometry In the left panel, the 
Kohn-Sham LDA band structure is compared with the 
scGW bands. Wc observe that scGW corrections upshift 
almost rigidly the lowest conduction states. Concern- 
ing the valence: (i) the dispersion of the S p-Cu d an- 
tibonding states at the top of the valence remains fairly 
unaltered, even if the overlap of LDA and quasiparticle 
wavefunctions shows remarkable variations in this region, 
(ii) The band width of the S p-Gu d bonding states (lo- 
cated between -3.5 and -6.7 eV) slightly increases, (hi) 
The bands associated to the (Zn,Sn)-S bond (between -8 
and -10 eV) are inverted and downshifted by about 2cV 
with respect to LDA. (iv) Also the Ss states are moved 
down by about 2 eV. In the right panel of Fig. Q] we show 
the same bands as obtained from GGA+U and HSE06 
calculations. In this case HSE06 bands are remarkably 
similar to scGW bands. This is not particularly sur- 
prising as the Hartrce-Fock mixing of HSE06 is particu- 
larly suited for materials with gaps of about 1-2 eV i 23 i 24 
As expected, GGA+U shifts down the states with Cu d 
character, thereby opening the gap to a reasonable value. 
However, it is evident from the figure, that the overall de- 
scription of the band dispersions is quite inaccurate. 

It is by now known that the band gap in CIGS materi- 
als is extremely sensitive to structural distortions^ 1 - and 
the remarkable stability of the band gap found experi- 
mentally can only be explained by compensating effects 
induced by intrinsic defects In order to establish if a 
similar behavior is also found in CZTS compounds, we 
performed calculations for kesterite (see Fig. [2]) and stan- 
nite Cu2ZnSnS4 by varying the anion displacements. We 
also verified that sensible variations of a and c lattice pa- 



2 

> 



1 

A-^LDA+1.5 eV 




1 


• HSE06 








»-«scGW 








- 






A--""' 






u z =0.13 (exp 0.128) 


m' 








0.23 


0.235 


0.24 
u =u 

x y 


0.245 0.25 



FIG. 2. (Color online) Quasiparticle band gap vs. anion 
displacement u x — u y for kesterite Cu2ZnSnS4, using DFT- 
LDA (red triangles), HSE06 (green circles) and scGW (blue 
squares). DFT-LDA values are upshifted by 1.5 eV. The ver- 
tical (horizontal) shaded areas give the spread of experimental 
data for u (band gap). 
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FIG. 3. Band structures from scGW for (a) kesterite 
Cu2ZnSnS4, (b) stannite Cu2ZnSnS4, (c) kesterite 
Cu2ZnSnSe4, (d) stannite Cu2ZnSnSe4. 



ramcters have negligible effects on the gap, in analogy to 
the case of CIGS.— We observe a strong variation of the 
band gap with u, which is due to similar variations of 
both valence and conduction band edges. Note that the 
position of the top valence has important implications 
for the formation energies of charged defects. The slopes 
are substantially larger for scGW calculations than for 
DFT-LDA, and even than for HSE06. In fact, they are 
controlled by the screening, which is essential to include 
in a self-consistent way, as in scGW. 

In Fig. [3] we can see the scGW band structures of the 
four compounds of the CZTS family. For the kesterite 
Cii2ZnSnS4 and stannite Cu2ZnSnSe4 the experimental 
geometries were usedJ^ For the remaining compounds, 
the experimental anion displacements are not reported^ 
and it was therefore necessary to resort to the theoreti- 
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TABLE I. Band gaps and valence-band shifts with respect to 
LDA (in eV) for all the structures considered in this work. 





LDA GGA+U HSE scGW exp. 


CD F 

S3 AE V 


n oq n sk r w i fid. i d—i f 14 
0.00 -0.52 -0.81 -0.49 


1 o Eg 
AE V 


-0.24 0.34 0.94 1.02 0.8-1.0^ 
0.00 -0.27 -0.55 -0.31 


a) E s 


-0.11 0.63 1.27 1.33 1.4-1.5^ 
0.00 -0.53 -0.80 -0.42 


« Fi — 


-0.41 0.15 0.75 0.87 0.8-1.0±& 
0.0 -0.28 -0.56 -0.35 



cal structures. In view of the strong variation of the gap 
mentioned before, the choice of the theoretical frame- 
work for the geometry optimization is essential. In fact, 
we verified that LDA/GGA relaxed structures lead to 
unacceptably large errors (up to 40%) in the scGW band 
gap, essentially related to the error in u: Cu-(S,Se) and 
Zn-(S,Se) bond lengths are too small in DFT-LDA, while 
Sn-(S,Se) bond lengths are too large. The solution to this 
issue relies on the use of the HSE06 functional, that yields 
extremely accurate values for u and the cell parameters, 
thanks to a better description of the localized states par- 
ticipating in the bonds. The major difference among the 
four band structures of Fig. [3] is the width of the band 
gaps, which are in excellent agreement with experimental 
results. As in the CIGS family, the Se compounds have 
a smaller gap with respect to the S ones. Furthermore, 
stannites have consistently smaller gaps than kestcrites. 

Finally, in table U we show the band gaps and the top 
valence shifts AE V with respect to DFT-LDA, obtained 
by aligning the average electrostatic potential in the dif- 
ferent theoretical schemes. The validity of this align- 
ment scheme is discussed in Ref. [25|. Band edge-shifts 
are essential quantities to determine the position of de- 
fect levels in the gap and band offsets at interfaces. The 
most striking feature is that, even if HSE06 gives very 
good band gaps, the top valence and bottom conduction 
bands are systematically too low. These results suggest 
that the fact that hybrids (and more in general tuned 
hybrids) can give good gaps, does not mean that they al- 
ways are able to reproduce correctly band-edge shifts^ 
We observe also that the valence band shift in GGA+U is 
close to the one obtained using scGW, despite the under- 
estimation of the gap. Note however, that GGA+U does 
not account for the dependence on u of the valence-band 
shift, that is particularly strong in this family. 

In conclusion, we obtained from many-body calcula- 
tions a very accurate description of the electronic prop- 
erties of the CZTS family. We proved that the HSE06 
hybrid functional offers a good compromise between ac- 
curacy and computational cost, yielding good gaps and 
relaxed structures. Nevertheless, it cannot assure a re- 
liable description of the valence and conduction contri- 
butions to the band-edge corrections, which has decisive 
implications for the position of the defect levels in the 
corrected band gap. 



Furthermore we found a remarkable similarity between 
the electronic properties of the stannites and kesterite 
compounds, and to a large extent between these and the 
CIGS family. This once more points to the direction 
that CZTS materials are excellent candidates to replace 
the more costly CIGS in the absorbing layer of thin film 
photovoltaic cells. 
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